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The ubiquitous inositol 1,4,5-trisphosphate (InsPs) receptor (InsPsR) channel, localized primarily in the endoplas- 
mic reticulum (ER) membrane, releases Ca^^ into the cytoplasm upon binding InsPs, generating and modulating 
intracellular Ca^^ signals that regulate numerous physiological processes. Together with the number of channels 
activated and the open probability of the active channels, the size of the unitary Ca^^ current (ica) passing through 
an open InsPjR channel determines the amount of Ca^* released from the ER store, and thus the amplitude and 
the spatial and temporal nature of Ca^* signals generated in response to extracellular stimuli. Despite its signifi- 
cance, ica for InsPjR channels in physiological ionic conditions has not been directly measured, fiere, we report 
the first measurement of Zca through an InsPjR channel in its native membrane environment under physiological 
ionic conditions. Nuclear patch clamp electrophysiology with rapid perfusion solution exchanges was used to study 
the conductance properties of recombinant homotetrameric rat type 3 InsPsR channels. Within physiological 
ranges of free Ca^* concentrations in the ER lumen ([Ca^^jER), free cytoplasmic [Ca'^] ([Ca^^Ji), and symmetric 
free [Mg^^] ([Mg^^]f), the «ca-[Ca^^]ER relation was linear, with no detectable dependence on [Mg^^]f. z'ca was 
0.15 ± 0.01 pAfor a filled ER store with 500 \iM [Ca^^lsR. The fca-[Ca^^]ER relation suggests that Ca^^ released by an 
InsPsR channel raises [Ca'^Ji near the open channel to ^13-70 \\M, depending on [Ca^^lsR. These measurements 
have implications for the activities of nearby InsPj-liganded InsPgR channels, and they confirm that Ca^^ released 
by an open InsPjR channel is sufficient to activate neighboring channels at appropriate distances away, promoting 
Ca^Mnduced Ca^^ release. 



INTRODUCTION 

Modulating cytoplasmic free Ca^^ concentration ( [Ca^*] i) 
is a ubiquitous intracellular signaling pathway that regu- 
lates numerous cellular physiological processes, includ- 
ing apoptosis, gene expression, bioenergetics, secretion, 
immune responses, fertilization, muscle contraction, syn- 
aptic transmission, and learning and memory (Clapham, 
1995; Berridge et al., 2000; Bootman et al., 2001; Braet 
et al., 2004; Randriamampita and Trautmann, 2004; 
Cardenas et al., 2010b) . The inositol 1,4,5-trisphosphate 
(InsPs) receptor (InsPsR) , a transmembrane protein lo- 
calized mainly at the ER in all animal cell types, plays a 
central role in this [Ca^^Ji signaling pathway (Taylor and 
Richardson, 1991; Bezprozvanny and Ehrlich, 1995; 
Furuichi and Mikoshiba, 1995; Patterson et al., 2004; 
Foskett etal., 2007; Joseph and Hajnoczky, 2007). In re- 
sponse to extracellular stimuli, phosphatidylinositol 
4,5-bisphosphate in the plasma membrane is hydrolyzed 
to generate InsPs (Berridge, 1993). InsPj rapidly dif- 
fuses through the cytoplasm to bind to the InsPjR and 
activates it as an intracellular Ca^* channel to release 
Ca^"^ stored inside the lumen of the ER into the cyto- 
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plasm, generating diverse local and global [Ca ]i sig- 
nals (Berridge, 1997). 

The InsPsR-mediated Ca'* flux from the ER store in 
response to various extracellular stimuli is Na icn Po, 
where A^a is the number of InsPgR channels activated, id 
is the unitary calcium ion current passing through an 
individual open InsPgR channel, and Pq is the open 
probability of the active InsPjR channels. The amount 
of Ca^* released, and therefore the amplitude and spa- 
tial and temporal nature of the [Ca^*]i signals generated, 
is directly dependent on «ca (Berridge, 1997; Bootman 
et al., 1997). Furthermore, the Pq of InsPgR channels is 
regulated by [Ca^^Ji with a biphasic dependence: at low 
concentrations, Ca^* activates the channel and increases 
its Po, whereas at higher concentrations, Ca^* inhibits 
the channel (Foskett et al., 2007). Consequently, z'ca 
also affects indirectly the amount of Ca^* released by 
regulating the Po of the activated channel itself, as well 
as that of nearby surrounding channels. Therefore, 
the meaisurement of ic^ in ionic conditions similar to 
those that exist physiologically in cells is critical for the 
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understanding of the mechanisms regulating this im- 
portant signaling pathway. 

Although InsPjR channel activity level (Pq) and the 
number of channels activated (A^) under various physio- 
logical conditions have been studied previously by 
electrophysiological methods, especially single-channel 
nuclear patch clamp experiments in various configura- 
tions (Foskett et al., 2007) , the unitary Ca^^ current (?ca) 
passing through an open InsPgR channel has not been 
characterized, primarily as a consequence of technical 
difficulties. Here, we measured the z'ca of recombinant 
homotetrameric rat type 3 InsPgR channels under physio- 
logical ionic conditions and studied its single-channel 
conductance properties. 

MATERIALS AND METHODS 

Nucleus isolation and nuclear patch clamp electrophysiology 

The generation and maintenance of DT40-KO-r-InsP3R-3 cells 
(mutant cells derived from chicken B cells with the endogenous 
genes for all three InsPgR isoforms knocked out and then stably 
transfected to express recombinant rat type 3 InsPsR) were de- 
scribed in Mak et al. (2005). Nuclear patch clamp experiments 
were performed using nuclei isolated from DT40-KO-r-InsP3R-3 
cells as described previously (Mak et al., 200.5). Excised nuclear 
membrane patches in the luminal side-out (lum-out) or cytoplas- 
mic side-out (cyto-out) configuration were obtained from isolated 
nuclei (Mak et al., 2007) using protocols analogous to those used 
to obtain inside-out or outside-out excised patches in plasma mem- 
brane patch clamp experiments. The solution around the excised 
nuclear membrane patch was rapidly switched multiple times using 
a solution-switching setup described in Mak et al. (2007). 

InsPsR channel current traces were acquired at room tempera- 
ture as described previously (Mak et al., 1998), digitized at 5 kHz, 
and anti-aliasing filtered at 1 kHz. Data analysis and the fitting of 
channel current-voltage curves were performed using Igor-Pro 
software. All electrical potentials were measured relative to the 
bath electrode. 

Experimental solution composition 

In cyto-out experiments, pipette solutions contained 140 mM KCl, 
10 mM HEPES, pH to 7.3 with KOH, 0.5 mM NagATP, and 10 pM 
InsPs, with free [Ca^*] ([Ca^Tf) buffered to 3 pM by 0.5 mM 
5,5'-dibromo l,2-bis(o-aminophenoxy) ethane-Af,A^,A'',A''-tetraacetic 
acid (dibromo BAPTA) and 0.2 mM CaClg. Perfusion solutions on 
the cytoplasmic side of the channel contained 10 mM HEPES, pH 
to 7.3 with KOH, 0.5 mM NajATP, and either 140 niM KCl with no 
InsPs or 70 mM KCl with 10 pM InsPj, with [Ca'+Jr buffered to 
3 pM by 0.5 mM (2-hydroxyethyl) ethylenediaminetriacetate and 
0.22 mM CaClj (Mak et al., 2005). InsPj and Na,ATP were in- 
cluded in the pipette solution to confirm that the cyto-out config- 
uration was properly achieved. Because of the presence of InsPs 
and ATP in the pipette solution, observation of InsPjR channel 
activity before solution switching would be evidence that the lum- 
out configuration was erroneously obtained. 

The same pipette solution was used in lum-out experiments to 
determine ion permeability ratios of the InsPgR channel. The per- 
fusion solution used to determine Pk: ^ci contained 30 mM KCl, 
no mM NMDG chloride, and 10 mM HEPES, pH to 7.3 with 
NMDG, with [Ca^*]f buffered to 300 nM by 0.5 mM l,2-bis(o-ami- 
nophenoxy) ethane- A',A',W,A''-tetraacetic acid (BAPTA) and 0.2 mM 
CaClg. The perfusion solution used to determine Pk- ^Na con- 
tained 140 mM NaCl and 10 mM HEPES, pH to 7.3 with NaOH. 



The perfusion solution used to determine Px'. Pk (X = Mg"*, Ca^*, 
St^\ or Ba^*) contained 140 mM KCl, 10 mM HEPES, pH to 7.3 
with KOH, and 10 mM XClj. The perfusion solutions contained 
no Na2ATP or Ca^"^ chelator. No CaClj was added to the perfusion 
solutions used to determine permeability ratios of Mg^*, Sr^*, and 
Ba^* versus K*, so free Ca^* in those solutions were from contami- 
nants in the water and salts used to make the solutions. Based on 
the purity assays of the salts (Sigma-Aldrich) and induction-coupled 
plasma mass spectr-ometiy assay (Mayo Medical Laboratory) of de- 
ionized water samples, [Ca^*]f in those solutions were -^5-8 pM. 

In lum-out experiments to determine the ica through the 
InsPsR-S channel, all solutions contained 10 mM HEPES, pH to 
7.3 with KOH. Unless stated otherwise, all solutions used in the 
same experiment contained the same concentrations of MgClj 
(0, 0.5, or I mM), KCl (140 or 40 mM), and potassium methane- 
sulfonate (KCH3SO3; 0 or 100 mM). Pipette solutions also con- 
tained 0.5 mM Na2ATP and 2 or 10 pM In.sP3 to activate the 
channel to P„ of ~0.5 (Mak et al., 2001a,b; Vais et al., 2010) , with 
[Ca^*]f buffered to either 70 nM by 0.5 mM BAPTA and 0.06 mM 
CaCls, or to 3 pM by 0.5 mM dibromo BAPTA and 0.2 mM CaClj. 
The same solution (without NajATP or InsP3) was used as perfu- 
sion solution for symmetric ionic conditions. Although [Na*] 
was not symmetric because of the presence of 0.5 mM NajATP 
only in the pipette solutions, the extra Na* is <I% of the amount 
of monovalent cations, and its contribution to channel current 
is negligible. Perfusion solutions used for asymmetric ionic con- 
ditions contained 0.3, I, or 2 mM CaClj, with no Ca^* chelator, 
Na,ATP, or InsPj. 

[Ca^lf of <100 pM (buffered by various Ca^* chelators) was 
confirmed by fluorimetry. [Ca^'^Jf in solutions without Ca^* chela- 
tors was calculated using activity coefficients (see next section). 

Evaluating ion permeability ratios 

According to the general Goldman-Hodgkin-Katz current equa- 
tion (Lewis, 1979), the permeability ratio iY/PsforY, one of the 
charge-carrying permeant ion species, can be evaluated from the 
reversal potential ( V^ev) of the channel as: 



Pv:P. 



[Y],-[Y]„exp(-Zyi'X„/i?r) 
l-exp(-ZYpK.„//?r) 

2 [x]i-[x]„exp(-z^pv;„/pr) ] 

l-exp(-Zxn;„/PT) 



(1) 



where X represents other permeant ion species present; Px is the 
permeability of X through the channel; Zx and zy are the valence 
of X and Y, respectively; [X];, [X]„, [Y];, and [Y]o are the activities 
of X and Yin the pipette and bath solutions, respectively; Pand R 
are the Faraday and gas constants, respectively; and Pis the abso- 
lute temperature, provided that all of the Px/Pk values are known. 
The activities of all ion species were calculated from activity co- 
efficients, except for [Ca^*] < 100 pM, which was determined 
by fluorimetry. 

The activity coefficients of KCl and NaCl were calculated using 
the Debye-Huckel equation based on data reported by Hamer 
and Wu (1972). The activity coefficients of CaClg in a I40-mM 
KCl solution (0.546-0.534 for 0.3-10 mM CaClz) were derived 
by interpolation using data reported by Butler (1968), assuming 
that CaClj activity coefficients are similar in 140-mM KCl and 
NaCl solutions. The activity coefficients of CaClj were used for 
other alkaline earth metal halides (MgCl2, BaClj, and SrCl2) in 
solutions with the same concentrations because activity coeffi- 
cients of the alkaline earth metal halides differ by <5% in con- 
centrations when activity coefficients are «0.54 (Goldberg and 
Nuttall, 1978). 
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RESULTS 

Conductance properties of homotetrameric reconnbinant 
rat InsPsR-S channels 

Although it has been a well-established approach to 
study the single-channel properties of the InsPgR chan- 
nel in its native membrane environment by performing 
patch clamp electrophysiology on isolated nuclei (Foskett 
et al., 2007) , the molecular composition of the channels 
observed in many of those studies could not be defini- 
tively ascertained. This is because the channels studied 
were either the endogenous channels (Mak and Foskett, 
1994, 1997; Marchenko etal., 2005; lonescu etal., 2006) 
with the possibility of different channel isoforms expressed 
and with possible alternative splicing (Foskett et al., 2007) , 
or recombinant channels expressed in cells with a non- 
zero level of endogenous InsPjR expression (Mak et al., 
2000; Boehning et al., 2001a), so that heteroligomeric 
channels of recombinant and endogenous InsPsRs could 
have been formed (Joseph et al., 1995; Mak et al., 2000) . 

To avoid possible variability that such heterogeneity 
might introduce, we ensured that only homotetrameric 
InsPsR channels with known identical amino acid se- 
quences were studied by using a stably transfected cell 
line (DT40-KO-r-InsP3R-3) derived from mutant DT40- 
InsPsR-KO cells (Sugawara etal., 1997) that have all en- 
dogenous genes for the three InsPgR isoforms knocked 
out, with only recombinant rat type 3 InsPsR (InsPsR-S) 
expressed (Mak et al., 2005; Li et al., 2007) . During our 



extensive experience working with this cell line (Mak 
et al., 2005; Foskett and Mak, 2010) , and in hundreds of 
nuclear membrane patches (in both cyto-out and lum- 
out configurations), only one kind of channel with a 
conductance >100 pS in symmetric 140 mM KCl was de- 
tected. The identity of these channels as recombinant 
rat InsPsR-S channels was confirmed by their sensitivity 
to activation by cytoplasmic InsPs (Fig. 1 A). 

Averaged over reasonable intervals (>1 s), the rat 
InsPjR-S channel dwells in a main open conductance 
state >95% of the time it is open, whereas it occasionally 
exhibits brief substates of lower conductances (Fig. 1, 
A and B). The substates are not a result of non-InsPsR 
channels because in current records showing only one 
active InsPjR channel gating, the channel current level 
dropped to substate levels from the main open state di- 
rectly without channel closing (Fig. 1 B, arrowheads) . 
Substates have also been observed for other InsPsR 
channels (endogenous or recombinant) in other cell 
systems (Watras et al., 1991 ; Mak and Foskett, 1997; Mak 
et al., 2000; Boehning et al., 2001a; lonescu et al., 2006) . 
To determine the conductance of the main open state 
of the channel in symmetric 140-mM KCl solutions 
([K"]f = [Cr]f 104 mM) in the absence of Mg'^, currents 
through excised lum-out membrane patches (/) were 
recorded as the applied potential (T4pp) was ramped 
(Fig. 1 C). Slope conductance of the channel (gch), 
evaluated as the difference between the slopes of the 
fits to open- and closed-channel current (/open and /dosed, 
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Figure 1. Current records of homotetrameric recombinant type 3 InsPgR channels in membrane patches excised from nuclei isolated 
from DT40-KO-r-InsP3R-3 cells. (A) Activation by InsPs of channels in a cyto-out nuclear membrane patch. Concentrations of InsPs and 
free Ca^* in the perfusion solution on the cytoplasmic side of the channel are indicated by color bars at top. \ipp = 30 mV. Channels 
opened only in the presence of cytoplasmic InsPj. Gray arrowheads indicate some of the occasions when one of the active channels 
entered a conductance substate. In this and subsequent current traces, black arrows indicate the membrane current levels when all 
InsPjR channels in the membrane patch were closed. (B) Single InsPgR-S channel current records observed in different lum-out nuclear 
patches in symmetric 140 mM KCl and 3 pM [Ca^^][. Pipette solution contained 10 pM InsPg. V^p^ = 40 mV. Gray arrowheads indicate 
when the channels entered different conductance substates with different conductance values. (C) Representative /-14pp plot for a single 
InsPjR-S channel in the same experimental conditions as in B. Kpp was ramped from —25 to 15 mV in 1 s. 20 ramps were analyzed. Data 
points selected for I^p^.^ and /dosed are plotted in green and pink, respectively. Green and red lines represent linear fits to /^pen and /dosed 
data points, respectively. 
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respectively) data, was 545 ± 7 pS (w = 16) . This and all 
subsequent open-channel current measurements were 
not affected by the presence of substates because atypi- 
cal current data arising from them were excluded when 
data points were selected for /-14pp fits. 

To properly design an experimental approach to 
measure the Zca passing through a single open InsPgR-S 
channel under physiological ionic conditions, it wais use- 
ful to obtain a maximum estimate of the size of the cur- 
rent using the Goldman-Hodgkin-Katz current equation 
(Hille, 2001). To obtain such an estimate requires knowl- 
edge of the permeability values of all permeant ionic 
species present. To evaluate the permeabilities of vari- 
ous ions through the InsPsR-S channel, /-T^pp data were 
recorded as T4pp was ramped during a series of lum-out 
patch clamp experiments in which the excised mem- 
brane patches were exposed to asymmetric ionic con- 
ditions (Fig. 2). After correcting for liquid junction 
potentials (Neher, 1995), the reversal potential (Kev) of 
the channel was determined as the 14pp at the intersec- 
tion of the /open-^4pp and /dosed- ''^pp fits. Permeability ra- 
tios for various ions were calculated from T^ev using the 
general Goldman-Hodgkin-Katz current equation (Eq. 1) . 
The permeability ratio for versus CP (Pk-Pca) was de- 
rived using asymmetric ionic solutions containing only 
two permeant ionic species: and CI (Fig. 2 A) . That 
value was then used to calculate the permeability ratios 
of other cations (Mg^^, Ca^^, as well as Ba^^, Sr'^, and 
Na^) versus from Kev measured in asymmetric ionic 
solutions containing three permeant ionic species 
(Fig. 2, B-F) , assuming that Pk:Pci is the same in all 
ionic conditions used. These measurements indicated 

that Pca-/sr-/Ba-/Mg-/Na-/K-/ci = (15.2 ± 0.6):(13.2 ± 

0.7):(11.8 ± 0.5):(10.2 ± 0.3):(1.24 ± 0.003):1:(0.27 ± 
0.01) (n = 3 for each ratio). This same InsPsR perme- 
ability ratio sequence wais also observed in other nuclear 
patch clamp experiments (Mak and Foskett, 1998; Mak 
et al., 2000; Boehning et al., 2001a; lonescu et al., 2006) . 

/'ca through InsPsR-S under physiological conditions 
Although the InsPsR channel has a large gch in solutions 
with physiological KCl concentrations, the current mea- 
sured under those conditions is mostly carried by 
ions driven through the channel by T{,pp. Because of the 
absence of a significant voltage across the ER membrane 
(Beeler et al., 1981; Marhl et al., 1997), Ca'* ions are 
driven through open InsPjR channels under physiologi- 
cal conditions by the difference between [Ca^*]f in the 
ERlumen ([Ca^^lER) and that in the cytoplasm ([Ca^^Ji). 
[Ca'^lER observed in various cell types is approximately 
hundreds of micromolars (Bygrave and Benedetti, 1996; 
Yu and Hinkle, 2000; Palmer et al., 2004), so there are 
significantly fewer Ca'* ions than K* ions to carry the 
current. A simple calculation based on the Goldman- 
Hodgkin-Katz current equation suggests that, in the ab- 
sence of transmembrane voltage, with [Ca"*]ER = 1 mM 



and [Ca^*]i = 70 nM, ica is '^3 pA for a channel with g-h = 
545 pS in 140 mM KCl and Pc^-Pk-Pci = 15:1:0.27. How- 
ever, the actual z'ca is expected to be substantially smaller 
because divalent cations act as permeant blockers that 
reduce gch (Mak and Foskett, 1998; Mak et al., 2000). 
Despite their higher permeabilities, divalent cations 
bind strongly to site (s) in the channel pore so that they 
pass through the channel significandy more slowly than 
monovalent cations. In addition, the high [K*] in the 
cytoplasm and ERlumen (^140 mM) and the relatively 
weak selectivity of the InsPsR channel for Ca^* over K* 
{Pch'-Pk = 15) make K* a potentially effective competing 
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Figure 2. Representative /-Kpp plots for single InsPgR-S channels 
in excised lum-out nuclear membrane patches under asymmetric 
ionic conditions used to determine channel permeability ratios 
for various ions. Compositions of pipette and perfusion solutions 
are described in Materials and methods. Selected data points for 
/open and /closed are plotted in green and pink, respectively. Red 
lines are linear fits to /dosed" ^pp data points, and green ones are 
polynomial fits to /open-Kipp data points (linear for B; quadratic for 
A, E, and F; cubic for C and D) . VJev at the intersection of the /open- 
IJpp and /closed" Kpp fits are marked by dashed lines and tabulated. 

(A) Eight 14pp ramps in one of four experiments to determine 
/k:/ci- Junction potential corrected (Inunction) = —9.3 ± 0.4 mV. 

(B) 15 14pp ramps in one of three experiments to determine 
Pn^-Pk- ^juncdon = "l-S ± 0.2 iiiV. (C) 10 K,pp ramps in one of 
three experiments to determine Pug'-pK- Ijuncdon = —0.7 ± 0.4 mV. 
(D) 43 Kpp ramps in one of three experiments to determine 
P^-Pk- ^unenon = -1-4 ± 0.2 mV. (E) 35 ramps in one of 
three experiments to determine Pb^.Pk. Vjuncnon = —2.6 ± 0.3 mV. 
(F) 16 14pp ramps in one of three experiments to determine Psr^^K- 
Kuncdon = -2.2 ± 0.8 mV. 
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ion to also reduce ?ca- «ca might be further reduced as a 
result of competition from free Mg^^, which exists at sig- 
nificant levels (^200 pM-1.1 mM) in both the cytoplaism 
and ER lumen (Halvorson et al., 1992; Morelle et al., 
1994a,b; Silverman et al., 1994; Golding and Golding, 
1995; Singh and Wisdom, 1995; Tashiro and Konishi, 
1997), and has a permeability similar to that of Ca^* 
(Pca:PMg= 1.5). 

To measure z'ca, we monitored the transmembrane 
current / through excised lum-out nuclear membrane 
patches containing a single active InsPgR channel. The 
capability to switch the perfusion solutions around the 
excised patch was used to minimize systematic measure- 
ment errors by calibrating the patch clamp setup for 
each individual experiment. The lum-out membrane 
patch was first exposed to a perfusion solution with the 
same ionic composition as the pipette solution (sym- 
metric ionic conditions) . /-T4pp data were recorded with 
l{pp ramped from —20 to 10 mV, and data points for 
^pen-Kpp and /closed- Kipp wcrc fitted linearly (Fig. 3 A). 
With no permeant ion concentration gradient across 
the excised nuclear membrane patch, the current pass- 
ing through the InsPaR channel (/open-^ciosed) and the 
leak current across the nuclear membrane patch (/dosed) 
must both be zero when T4pp = 0. Thus, the intersection 
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Figure 3. Measuring ica in physiological ionic conditions. 
(A) Linear fits (green and red lines) to selected /open and /dosed 
data (plotted in green and pink, respectively) from 25 l^pp ramps 
recorded from an excised lum-out nuclear membrane patch con- 
taining one active InsPsR channel under symmetric ionic condi- 
tions, with pipette and perfusion solutions containing 140 mM 
KCl, 0.5 mM MgCl2, and 3 pM [Ca^*]f. Pipette solution contained 
2 }iM InsPa. (B) Graph of the fitted /open-^pp and /dosed- ^app lines 
in the /-14pp region marked by the black rectangle in A. Zero-cur- 
rent level (black dotted line) established at the intersection of the 
^-14pp fits at I4pp = 0 (marked by orange arrowhead). (C) Linear 
fits to /open and /dosed data (same convention as in A) from 25 ^^pp 
ramps recorded for the same membrane patch under asymmetric 
ionic conditions, with perfusion solution containing 2 mM CaCl2. 
/and \4pp ranges in A and C are the same. (D) Graph of the fitted 
/-Kpp lines in the same /-14pp region with the same zero-current 
level as in B. \4pp = 0 (marked by orange arrowhead) at the inter- 
section of the /dosed- Kipp line and zero-current level, ica is /open at 
14pp = 0 (marked by blue arrow) . 



of the /open-^4pp and /dosed- ''^pp lines established the zero- 
current level for the patch clamp current recording sys- 
tem during each experiment (Fig. 3 B). 

The perfusion solution was then switched to one with 
higher [Ca^*]f (asymmetric ionic conditions). Higher 
[Ca'^] ER reduced both the InsPgR g^i, as a result of an in- 
crease in permeant divalent cation block, and the leak 
conductance (Fig. 3 C). However, the zero-current level 
of the recording system was not affected by the perfu- 
sion solution switch. For the leak current (/dosed) through 
the excised membrane patch, T^ev = 0 because there are 
no significant concentration differences for the major 
ionic components in the pipette and perfusion solu- 
tions, and CI , which have significantly higher mobil- 
ity across the membrane than divalent cations (Beeler 
et al., 1981; Marhl et al., 1997). Therefore, l^pp = 0 at 
the point where /dosed = 0 (Fig. 3 D). The open-channel 
current (/open) at V^pp = 0 is driven only by the [Ca^^Jf 
gradient across the channel ( A[Ca^*] = [Ca^^ljn -[Ca^*]; ) 
and therefore is Zca (Fig. 3 D) . 

We first measured Zca in symmetric 140 mM KCl with 
no Mg^^. To cover the range of reported [Ca^^lER (Bygrave 
and Benedetti, 1996; Yu and Hinkle, 2000; Palmer et al., 
2004) , the luminal side of the excised membrane patch 
wais exposed to perfusion solutions with [Ca^^Jf = 160 pM, 
550 pM, and 1.1 mM. [Ca^^Ji in the pipette solution was 
fixed at 3 pM to keep channel Pq high. The observed 
magnitude of ica was well described by a linear relation 
with A[Ca"*] (Fig. 4), as predicted by the Goldman- 
Hodgkin-Katz current equation. However, the value for 
Pca (1.5 X lO"^'^ m^s"^) derived from the slope of the z'ca 
versus A[Ca^^] line is an order of magnitude smaller than 
that (1.5 X 10"" m'^s"^) estimated using g-d, in 140 mM 
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Figure 4. zca's through the InsPjR channel under various asym- 
metric ionic conditions. Error bars show SEM. The number of 
measurements performed for each ionic condition is tabulated 
next to the corresponding data point. The line is the linear fit to 
the data points in 140 mM KCl, 0 mM [Mg^Tj, and .3 pM [Ca^*];, 
with slope of 0.30 ± 0.02 pA/mM. 
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KCl and the measured permeability ratios of Pca^^k^^ci- 
This discrepancy results from interactions between the 
channel and the permeant ions (K^ and Ca^^) and among 
permeant ions in the channel that are not considered 
in the Goldman-Hodgkin-Katz equation. 

To assess the effect on ica of physiological concentra- 
tions of free Mg^^, ranging from 200 \iM to 1.1 mM in 
both the cytoplcism and ER lumen (Halvorson et al., 1992; 
Morelle et al., 1994a,b; Silverman et al., 1994; Golding 
and Golding, 1995; Singh and Wisdom, 1995; Tashiro 
and Konishi, 1997), we measured «ca in symmetric 270 
or 550 pM [Mg2+]f with [Ca^1i = 3 pM and [Ca^+lER = 
550 pM or 1.1 mM. Interestingly, «ca observed in the 
presence of physiological [Mg^^]f was not significandy 
different from Zca measured in 0 [Mg~*]f (Fig. 4). 

Under normal physiological conditions, [Ca^^jER is 
significantly higher than [Ca^^] i, so Zca should have little 
dependence on [Ca^^Jj. Indeed, we observed no statisti- 
cal difference between zci, for [Ca^^Ji at resting (70 nM) 
and activating (3 pM) levels (Fig. 4). 

Free chloride ion concentration [CI ]f is significandy 
lower than [K^]f in the cytoplasm as a result of the Gibbs- 
Donnan effect of negative charges in cytoplasmic proteins 
(Foskett, 1990). To verify if physiological cytoplasmic 
[CI ]f significantly affects ica, we measured Jca with pi- 
pette and perfusion solutions containing symmetric 
100 mM potassium methanesulfonate (KCH3SO3) , 40 mM 
KCl ([K"]f = 104 mM, [Cr]f ^ 30 mM), and 270 pM 
[Mg^lf, with [Ca^'li = 3 pM and [Ca'+JER = 550 pM or 
1.1 mM. Again, ic-^ observed wais not significandy different 
from that observed in 140 mM KCl (Fig. 4), indicating 
that zca is to a large extent independent of [CI ]f under 
physiological ionic conditions. Interestingly, g-h = 293 ± 4 
pS in symmetric 270 pM [Mg^^+Jf and 30 mM [CHf, the 
same as that observed in symmetric 270 pM [Mg^*]f and 
104 mM [Cr]f (292 ± 5 pS). This correspondence proba- 
bly arises because decreasing [CI ]f reduces CI current 
through the channel and also reduces CI competition 
with K* to move through the channel, with the two op- 
posing effects on InsPsR gch cancelling each other out. 

These results indicate that, in physiological ionic con- 
ditions, with symmetric 104 mM [K^]f, 30-104 mM 
[Cr]f, 0-550 pM [Mg'+]f, 70 nM to 3 pM [Ca^^];, the Zca 
for rat homotetrameric type 3 InsPsR channel in native 
ER membrane is (0.30 ± 0.02 pAmM^^) x [Ca^^ER- This 
is equivalent to having an asymptotic Ca'^ slope conduc- 
tance {gch as V^pp in the ER relative to the cytoplasm— >oo) 
of (23.4 ± 1.6 pS mM"^) x [Ca'^jER. For a filled store 
with 500 pM [Ca^^lER, «ca = 0.15 ± 0.01 pA and the asymp- 
totic Ca^^ slope conductance is 11.7 ± 0.8 pS. 

Block of InsPsR channel conduction by permeant 
divalent cations 

Another indication of the inadequacy of the Goldman- 
Hodgkin-Katz equation to describe ionic flow through 
the InsPsR channel is the substantial reduction of gch by 



Ca'^ or Mg'* (in millimolar concentrations) on the lumi- 
nal or cytoplasmic side of the channel, despite their high 
permeability ratios (Fig. 5 A) . Similar partial reduction of 
gch by permeant divalent cations (Mg^*, Ca'^, and Ba^*) 
Wcis observed in endogenous type 1 InsFjR channels in 
Xenopus laevis oocyte nuclear membrane (Mak and Foskett, 
1998) and recombinant type 1 InsPsR channels in the 
plasma membrane of DT40-KO-r-InsP3R-l cells (Dellis 
et al., 2006). This reduction suggests that movement of 
divalent cations through the InsPgR channel is substan- 
tially slowed down by strong interactions between the 
ions and the pore, causing temporary block of ion flow 
through the channel (Hille, 2001). The channel conduc- 
tance gch at the reversal potential l^ev was systematically 
evaluated in experiments performed to measure zca in 
the presence of 3 pM, 160 pM, 550 pM, and 1.1 mM 
[Ca'^ER, and symmetric 0, 270 pM and 550 pM [Mg^+]f. 
Channel blocking effects of Ca^* and Mg'^ were not mu- 
tually exclusive, as the addition of Mg^^ further reduced 
gch already suppressed by Ca'*, and vice versa (Fig. 5 A) . 
The observed reduction of gch by Ca"* could be described 
by an empirical saturating partial inhibition equation: 

(2) 

where go is the channel slope conductance at 0 [Ca'^JsR, 
g»caER is the channel slope conductance at saturating 
[Ca^^] ER, and -KcaER is the half-maximal blocking [Ca er. 
go, g»caER, and -KcaER have different values in different 
[Mg^^]f (corresponding to the different curves in Fig. 5 A) . 
In the absence of Mg~^ (Fig- 5 red curve), Kcser = 
210 ± 20 pM, g»caER = 82 ± 8 pS, and go is 545 pS, the 
channel slope conductance with no divalent cations. 

To compare the capacities of Mg'* and Ca^* to block 
the channel, we measured gch with various [Mg^*] er, in 
0 [Mg^^^li and symmetric 3 pM [Ca^^Jt (Fig. 5 B, purple 
open squares), gch data could be fitted using an equa- 
tion similar to Eq. 2: 

(3) 

gci, ([Mg'^lER) = (,?„MgER[Mg'*]ER + ,?oJ*^MgER)/(^MgER + I-'^S^'^^r)- 

Eq. 3 describes the purple curve in Fig. 5, with go again 
being the channel conductance in the absence of diva- 
lent cations (545 pS). Interestingly, g»MgER = 81 ± 5 pS is 
very similar to g>=caER, suggesting that Ca^* and Mg'^ are 
equally efficacious in blocking the InsPsR channel. Fur- 
thermore, XivigER = 410 ± 20 pM, which is = 2 x i^^^ER, sug- 
gesting that the affinity for Ca^* of the site in the channel 
responsible for blockage by luminal divalent cations is 
twice that for Mg^^; i.e., Mg'* is half as potent in block- 
ing the InsPjR channel as Ca^^. 

Although inhibition of InsPsR gating by high [Ca^^Ji 
(>20 pM) (Mak et al., 2001b) prevents comparison of 
the effectiveness of Ca^* from the cytoplaismic and Itiminal 
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Figure 5. Slope conductance of an InsPaR-S channel at IJev in sym- 
metric 140-mM KCl solutions with various divalent cation concen- 
ti'ations. SEM for all data points are shown as error bars, some of 
which are smaller than the size of the symbols. The number of mea- 
surements performed for each ionic condition is tabulated next to 
the corresponding data point in the same color. The same symbols 
in the same color are used for the same set of data points plotted in 
the different graphs. (A) g.,, versus [Ca^*]ER in the presence of vari- 
ous symmetric [Mg^*]f as tabulated. Curves are fits to gih data by the 
empirical inhibition equation (Eq. 2), with different Xcser for vari- 
ous [Mg^*][. The red curve is for g-i, in 0 [Mg^^]f. (B) g^i, versus con- 
centrations of divalent cations in the cytoplasmic or luminal side of 
the channel as tabulated. The purple curve is described by Eq. 3, and 
the red one is described by Eq. 2. (C) versus equivalent divalent 
cation concentrations [X^*]eq= [Ca^*]ER-^ 0.5 x ([Mg^*]i4- [Mg^*]ER) 
for various combinations of [Ca^*]ER, [Mg^*]ER, and [Mg^*]i as tabu- 
lated. The red curve is described by Eq. 5. 



g-ch observed in 0.55 mM [Mg'*]i and 0 [Mg'^jER (Fig. 5 B, 
magenta filled square) was the same as that in 0 [Mg 
and 0.55 mM [Mg^^jER. This indicates that Mg^* blocks 
the InsPgR channel with equal potency firom either side. 
Thus, the site responsible for permeant divalent cation 
block is probably located inside the pore along the ion 
permeation pathway, with cations fi"om either side of the 
channel having similar access to the site. 

The data for g-h in symmetric [Mg'*]f (Fig. 5 B, black 
crosses) fall on the curve described by Eq. 2 (Fig. 5 B, 
red curve), so that 



^ch ([Mg'"],,^) = (g.CaER[Mg'"],^ + ff„^:CaER)/(^CaER + [Mg'*],,™). 



Using g-„caER = g»MgER = g=<., and Kc-^er = 0.5 KugEs., 



MgER 



^MgER+{[Mg'']ER+[Mg'"]J (4) 



suggesting that contributions to channel block by Mg^* 
on either side can be combined by simply summing the 
[Mg^*]f on the two sides as if all the Mg"* was on one 
side. This relation can even be extended to include 
[Ca^^lER- gch data observed in various combinations of 
[Ca^;]ER, [Mg^^];, and [Mg^+JsR (whether [Mg^^]; = 
[Mg^^lER or not), when plotted against the equivalent 
divalent ion concentrations [X^^Jeq defined as 
|[Ca2-^]ER + 0.5x([Mg2+]ER + [Mg2ni)}, are all well fitted bya 
curve similar to that described by Eq. 2 (Fig. 5 C, red 
curve); i.e., 

gch([Ca^"]ER.[Mg'"]ER,[Mg2"]i) 



[Ca2"]ER + 



0.5 X ([Mg^-^ ]er + [Mg'* I )} + goXc,ER 

^caER +{[Ca^"]ER + 0.5x([Mg2-]ER + [Mg^^],)} 

-'^'caER + [^^*]en 



(5) 



This result strongly indicates that channel block by dif- 
ferent permeant divalent cations is caused by binding to 
a unique, saturable site in the ion permeation pathway 
that is equally accessible from either side of the channel, 
with an affinity for Ca"* twice that for Mg^*. This site is 
probably located at or near the selectivity filter of the 
channel where the channel pore size is most restricted. 



DISCUSSION 



sides to block permeation, the effectiveness of Mg to 
block permeation from either side of the channel can 
be compared because physiological [Mg^^Ji has no 
significant effect on channel gating (Mak et al., 1999). 



Magnitude of /'ca through the insPsR channel 
in physiological ionic conditions 

In this study, the conduction properties of a homotetra- 
meric recombinant rat InsPgR-S channel were examined 
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under various ionic conditions, especially under physio- 
logical monovalent and divalent ion concentrations. 
Under physiological ranges of [K^]f, [CP]f, [Mg'*]f, 
[Ca'^li, and [Ca^^lER, «ca through an open InsPsR chan- 
nel depends only on [Ca^^lsR with a linear relation, = 
0.15 ± 0.01 pAfor a filled ER store with 500 pM [Ca^^lER- 
This value is compatible with the magnitude of Ca^^ flux 
estimated by the imaging of Ca^^ release events (puffs) 
of 0.4-2.5 pA (Sun etal., 1998) and 0.12-0.95 pA (Bruno 
et al., 2010) in Xenopus oocytes, where multiple active 
InsPsR channels were involved in generating a puff. The 
value is also in reasonable agreement with the values for 
«ca estimated from Ca^^ release events (blips and puffs) : 
-^0.4 pA (Shuai et al., 2006) and ^0.1 pA (Bruno et al., 
2010) for endogenous InsPsR-l in Xenopus oocytes, and 
^0.05 pA estimated for endogenous InsPsR-l chan- 
nels in human neuroblastoma SH-SY5Y cells (Smith and 
Parker, 2009) . Accordingly, our experimental approach, 
with the advantages of observing Ca^^ currents with sin- 
gle-channel resolution and rigorous control of ionic 
conditions on both sides of the channel, appears to ac- 
curately reflect the physiological behavior of InsPsR 
channels in intact cells. Furthermore, our measured 
value of ica is comparable to the z'ca of 0.1 pA (Swillens 
et al., 1999), 0.07 pA (Thul and Falcke, 2004), and 0.2 
pA (Shuai et al., 2008) , assumed in various efforts to nu- 
merically simulate Ca^^ release through InsPjR chan- 
nels, thus providing experimental support for the validity 
of those modeling efforts. 

icn for InsPjR measured here is comparable to but 
smaller than those determined for RYR channels, the 
other major family of intracellular Ca'^ release chan- 
nels, zca driven by 500 \iM [Ca'*]ER through purified 
amphibian type 1 RYR (RYRl) and mammalian type 2 
RYR (RYR2) channels reconstituted in artificial lipid 
bilayers was estimated to be 0.26 and 0.27 pA, respec- 
tively, in the presence of symmetric 150 mM KCl and 
1 mM MgCla (Kettlun et al., 2003). Under similar 
ionic conditions (140 mM KCl and 1 mM MgCl2), ?ca 
through an open rat InsPjR-S channel is 0.15 pA. 
However, as discussed in more detail below, if the lipid 
environment impinges on the conductance properties 
of the release channels, a comparison of the icn may 
not yet be possible because ?ca measurements for the 
two channel types were obtained in different mem- 
brane environments. 

InsPaR channel conductance 

gch of various InsPjR channel isoforms in symmetric 140 
mM KCl has been observed, mainly by nuclear patch 
clamp experiments (Mak and Foskett, 1998; Mak et al., 
2000; Boehning et al., 2001a; lonescu et al., 2006; 
Betzenhauser et al., 2009), to be comparable to but 
smaller than those for RYR channels reconstituted into 
lipid bilayers under the same KCl concentration: -750 
pS for RYRl (Wang et al., 2005) and 740 pS for RYR2 



(Lindsay etal., 1991). Comparisons of primary sequences 
of different InsPsR isoforms from various species with 
those of RYR isoforms reveal that a conserved se- 
quence GGGXGDX (amino acid residues 2545-2551 in 
rat InsPsR-l Sr Sir isoform [Mignery et al., 1990] and 
2472-2478 in rat InsPjR-S [Blondel et al., 1993]; X 
stands for I or V) between putative transmembrane he- 
lices 5 and 6 in the InsPjR pore-forming domain is 
highly homologous to the sequence GGGIGDE (amino 
acid residues 4895-4901 in human RYRl [Fujii et al., 
1991] and 4824-4830 in human RYR2 [Zorzato et al., 
1990]), which is conserved in all three RYR isoforms. 
Site-directed mutagenesis in the homologous RYR se- 
quence alters the conductance properties of RYR chan- 
nels (Zhao et al., 1999; Gao et al., 2000; Du et al., 2001; 
Chen et al., 2002; Wang et al., 2005), suggesting that 
those amino acids lie in or near the selectivity filter that 
determines, at least in part, the conductance of the RYR 
channel. Point mutations in the corresponding sequence 
in the rat InsPsR-l also altered its conductance proper- 
ties (Boehning et al., 2001b). A mutation changing 
the GGGVGDV sequence to GGG/GDV made it more 
similar to that of the RYR and increased InsPgR chan- 
nel conductance (Boehning etal., 2001b). Conversely, 
GGG/GDE to GGGVGDE (Gao et al., 2000) and GGGIGDE 
to GGGIGDQ (Wang et al., 2005) substitutions in the 
RYR sequence generated channels with reduced con- 
ductance. Furthermore, invertebrate InsPjR isoforms have 
a GGGIGDI sequence (Yoshikawa et al., 1992; Baylis 
et al., 1999; Iwasaki et al., 2002) that resembles the RYR 
sequence and have a higher single-channel conduc- 
tance (477 ± 3 pS; lonescu et al., 2006) than the verte- 
brate isoforms (360 - 390 pS; Mak and Foskett, 1998; 
Mak et al., 2000; Boehning et al., 2001a; Betzenhauser 
et al., 2009), which have a GGGVGDX sequence. Col- 
lectively, these observations strongly suggest that the 
GGGXGDX sequence in InsPjR is close to or forms part 
of the selectivity filter in the tetrameric channel, and is 
therefore a major factor that determines the conduc- 
tance properties of InsPjR channels. However, muta- 
genesis suggests that amino acids outside the GGGXGDX 
sequence also play a role in determining the conduc- 
tance of the InsPsR (unpublished data) and RYR (Gao 
et al., 2000; Du et al., 2001; Xu et al., 2006) channels, 
possibly through electrostatic effects to concentrate ions 
in the channel. 

Other factors also appear to contribute to the chan- 
nel conductance properties. Different single-channel 
conductances have been observed for the same recom- 
binant InsPsR-l isoform expressed in DT40-InsP3R-KO 
cells, depending on whether it was localized in the outer 
nuclear membrane (373 ± 2 pS; Betzenhauser et al., 
2009) or the plasma membrane (214 ± 17 pS; Dellis 
et al., 2006). Moreover, the conductance of the homo- 
tetrameric recombinant InsP3R-3 channel expressed in 
the outer nuclear envelope of DT40-KO-r-InsP3R-3 cells 
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measured here (545 ± 7 pS) is nearly 50% larger than 
that observed for the same recombinant InsPsR-S chan- 
nel expressed in the outer nuclear envelope of Xenopus 
oocytes (370 ± 8 pS; Mak et al., 2000). Conductance 
values ('^125 and 200 pS) smaller than the one we ob- 
served here were reported for the same channel in the 
same location from the same cell type (Taufiq-Ur-Rahman 
et al., 2009), but those smaller and variable conduc- 
tances were likely a result of heavy contamination by 
Mg2+ of the NagATP used (Rahman and Taylor, 2009) . 
The conductance of the r-InsPgR-S channels in the 
outer nuclear membrane of DT40-KO-r-InsP3R-3 cells is 
even larger than that for the invertebrate InsPjR (477 ± 
3 pS; lonescu et al., 2006), despite the more RYR-like 
putative selectivity filter sequence of the latter. Because 
all of these patch clamp studies of various InsPgR chan- 
nels were all performed in symmetric 140 mM KCl 
([K*]f = 104 mM), these observations suggest that be- 
sides the primary sequences of the InsPgR, other 
factor (s) — lipid environment, interacting proteins or 
peptides, etc., that are different in various expression 
systems (outer nuclear membrane vs. plasma mem- 
brane; DT40-KO-r-InsP3R-3 nucleus vs. Xenopus oocyte 
nucleus) — can affect the conductance of InsPsR chan- 
nels substantially. 

The factors that cause the differences of the InsP3R-3 
gch observed in the same rigorously controlled ionic 
conditions (symmetric 140 mM KCl) in different cellu- 
lar locations and different cells may also affect the size 
of Zca through InsP3R-3 channels in those contexts. 
Because ours is the first measurement of icn through an 
InsPsR channel under physiological ionic conditions, 
how «ca correlates with g^^ in various cell systems cannot 
be clearly gauged until similar measurements are made 
in different systems. It should be pointed out that al- 
though the ionic conditions ([Ca^^jER, [Mg^'^Jf, [K+], 
and [CI ]) under which zca was measured in this study 
were physiological or near-physiological, the ionic con- 
ditions under which g^h in symmetric 140 mM KCl were 
measured (in this study and all published reports) were 
nonphysiological because of the absence of divalent cat- 
ions, especially Mg'^, on either side of the channels. It is 
possible that the observed variability of gch is a result of 
the use of nonphysiological ionic conditions, as a study 
on the effects of divalent cations on InsPsR gVh suggested 
(Mak and Foskett, 1998) . This issue can be clarified with 
more direct measurements of ?ca for other InsPjR chan- 
nels in different cell systems. At this point, the applica- 
bility of the zca of homotetrameric type 3 InsPjR channels 
in DT40-KO-r-InsP3R-3 cells to other InsPjR channels 
and cell types should be considered judiciously with rec- 
ognition of its potential limits. However, this first direct 
measurement is of significant value for improving our 
understanding of the mechanisms of InsPsR-medi- 
ated Ca^* release and for modeling efforts to simulate 
Ca^^ signaling. 



Selective permeant ion block of the InsPsR channel 
by Mg2+ 

Physiological [Mg^^Jf (200 pM-l.l mM in both the cyto- 
plasm and ER lumen) substantially reduced gch of the 
InsPsR channel in symmetric 140-mM KCl solutions 
(Fig. 5) by acting as a permeant blocking cation. This 
indicates that under physiological conditions, Mg^^ sig- 
nificantly affects the passage of K^ through the InsPsR 
channel pore. Accordingly, it is not immediately obvi- 
ous why ica was not measurably different in the presence 
or absence of physiological [Mg~^]f (Fig. 4). Similarly, 
Zca's passing through RYR channels under physiological 
(-^500 pM) [Ca-^lER (Chen et al., 2003; Ketdun et al., 
2003; Gillespie and Fill, 2008) were not significantly af- 
fected by the presence or absence of 1 mM MgCl2. 

Because of the lack, to date, of a quantitative model 
to describe the conductance properties of the InsPsR 
channel pore, we attempt to explain qualitatively the 
apparently contradicting observations by using models 
developed to account for conductance properties of 
RYR channels, which have a putative selectivity filter se- 
quence, and therefore structure, highly homologous to 
those of InsPsR channels. 

In a barrier model that describes the RYR channel as 
a single-ion occupancy channel with four energy barri- 
ers in the selectivity filter (Tinker et al., 1992), the 
channel is blocked when any cation enters the vacant 
selectivity filter because the channel cannot accommo- 
date two cations simultaneously. Despite higher perme- 
ability of Mg~* than K*, because [K^]f in the cytoplasm 
and ER lumen is about two orders of magnitude higher 
than [Mg^^]f, K^ is likely to be the major permeant 
blocker of Zci, through the InsPsR channel. Further- 
more, because of the higher affinity of Ca^^ for the se- 
lectivity filter relative to that of Mg^^, as revealed by its 
more potent reduction of InsPsRg;!! (Fig. 5), the energy 
released as Ca^^ enters the channel selectivity filter from 
the luminal side can compensate for the energy re- 
quired to push an occupying Mg^* out the cytoplasmic 
side, especially if the Mg'* is bound to the side energy 
well close to the cytoplasmic end of the filter. Thus, it is 
possible that the id through the InsPjR channel is al- 
ready significandy suppressed by physiological [K^]f, 
such that Mg^^, with lower affinity for the selectivity 
filter than Ca^*, cannot reduce the magnitude of «ca fur- 
ther by a detectable amount. 

In Poisson-Nemst Planck models with (Gillespie et 
al., 2005) and without (Chen et al., 1997, 2003) density 
function theory, a cation-selective channel is described 
as one with a constricted selectivity filter region contain- 
ing negative charges from acidic amino side chains or 
carbonyl backbones. Cations move through the pore 
by electrodiffusion under electrical and concentration 
gradients, interacting with the channel both electro- 
statically and chemically. According to these models, 
concentrations of cations in the selectivity filter of the 
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RYR/InsPgR channel are higher than those in the bulk 
solutions because of the permanent negative charges 
from the Asp (in GGGXGDX in InsPgR) or Asp and Glu 
(GGGIGZ'A'in RYR) located in or near the selectivity fil- 
ter of the channel. This concentrating effect is stronger 
for Ca^^ and Mg^^ than for because of the more nega- 
tive chemical potentials of Ca'* and Mg^^ in the selectiv- 
ity filter as a result of either chemical interactions 
between the divalent cations and the channel (Chen et al., 
2003) , or because of the smaller ionic radii and there- 
fore smaller excluded volumes and higher charge 
densities of Ca^" and Mg'^ (Gillespie, 2008). The high 
divalent cation concentration in or near the selectivity 
filter in turn lowers [K^] there as a result of electro- 
static repulsion between the cations (Chen et al., 2003; 
Gillespie et al., 2005; Gillespie, 2008). This can cause 
the observed reduction of gch by physiological concen- 
trations of Mg^^ or Ca"*. In InsPgR channels, the con- 
served sequence GGGXGDX in or near the selectivity 
filter has fewer acidic residues than the corresponding 
sequence GGGIGDiifor the RYR channels. This may re- 
duce the permanent negative charge density in the 
InsPsR channel selectivity filter relative to that for the 
RYR channel and consequendy weaken the electrostatic 
interaction between the selectivity filter and the cations 
(Gillespie et al., 2005), making the effect of the more 
negative chemical potential of Ca'* than Mg^* in the fil- 
ter even more prominent in the InsPjR than RYR, as in- 
dicated by the substantial higher potency of Ca^^ than 
Mg^^ to reduce InsPjR g;h (Fig. 5). High [Ca^^Jf in the 
InsPsR selectivity filter will suppress entry of Mg'^ 
through electrostatic repulsion. With only physiological 
[Mg~*]f in the bulk solution around the channel, there 
may not be enough Mg^^ in the filter to significantly im- 
pede the flux of Ca^^. 

[Ca^^li at various distances from an open InsPsR channel 
as a result of ica 

The size of «ca through an open InsPjR channel in the 
ER membrane obviously directly affects the [Ca^^]; im- 
mediately surrounding it. Measurement of «ca in physio- 
logical ionic conditions now enables us to estimate the 
[Ca^*] i in the vicinity of an open InsPjR. For the estima- 
tion, the channel pore can be treated as a point source 
in a semi-infinite region (Smith, 1996). Within short 
distances from the open-channel pore (<15 nm), the 
concentrations and Ca^*-binding rates of endogenous 
cytoplasmic Ca^^ buffers are too low to significantly af- 
fect the local [Ca"^];. Thus, local [Ca^^]; in the immedi- 
ate vicinity of the open channel is mainly determined by 
the equilibrium between Ca^* flux through the channel 
and diffusion of unbound Ca^^ through the cytosol, and 
can be estimated with reasonable accuracy as 

[Ca^"'].(r)»[«c, /(zca^27tDr)] = l.l X [Ca^+J^^nm / r, (6) 



where r is the distance from the channel pore, D is the 
diffusion coefficient for Ca^* in the cytoplasm with no 
buffering, and Zca and /"are as defined in Eq. 1 (Smith, 
1996; Neher, 1998), using the value of ?ca measured 
here and D ^ 225 pm^s"^ (Allbritton et al., 1992). Ac- 
cording to cryo-electron microscopy (Jiang et al., 2002; 
da Fonseca et al., 2003; Hamada et al., 2003; Serysheva 
et al., 2003; Sato et al., 2004; Wolfram et al., 2010) and 
electron microscopy (Cardenas et al., 2010a) measure- 
ments, the radius of a tetrameric InsPjR channel in a 
plane perpendicular to the axis of the channel pore is 
^10-12 nm. There is no structural information con- 
cerning the locations of various Ca^^-binding sites in the 
channel relative to the pore, but it is reasonable to as- 
sume that the distance between the sites and the pore is 
<10-12 nm. Using 8 nm as an estimate of the distance 
between the channel pore and the activating and inhib- 
itory Ca^^-binding sites of the channel, within the physi- 
ological range of [Ca^+lER (100-500 pM), [Ca^^]; is 
sensed by the channel; i.e., [Ca'^Ji (8 nm) = 14-69 pM 
(Fig. 6) . This [Ca^*], level is reached within a short time 
(microsecond) after the opening of the channel (Neher, 
1998). With a filled ER store, [Ca'l; (8 nm) is high 
enough that even at saturating [InsPs] , Po of the open 
channel itself is significandy inhibited for most InsPsR 
isoforms that have been studied in single-channel ex- 
periments (Foskett et al., 2007). This can provide nega- 
tive feedback to terminate the Ca^* release. On the 
other hand, when the ER store is partially depleted, Ca^* 
released by the InsPsR channel may not be sufficient to 
raise [Ca'*]i high enough to suppress activity of the re- 
leasing channel to terminate the Ca^^ release. 

Farther away from the channel pore, [Ca^^] i is strongly 
affected by cytoplasmic Ca^* buffers. Estimates of con- 
centrations and Ca^^-binding rates of endogenous Ca^* 

[2+1 2+ 
Ca Jer Ca buffering 

500 |jM weak 




Distance from channel pore (nm) 

Figure 6. Spatial profiles of [Ca^*]; at various distances from an 
open InsPgR channel for different [Ca^*]ER with different cyto- 
plasmic Ca^*-buffering capacities. [Ca^*]i was calculated using 
Eq. 7 with characteristic length \ = 55 and 440 nm for strong and 
weak cytoplasmic Ca^"^ buffering, respectively. 
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buffer(s) (Bj and Ao,,, respectively) vary widely (Wagner 
and Keizer, 1994; Smith et al., 1996; Falcke, 2003; Shuai 
et al., 2008), so only a first-order estimate of [Ca^*]i(r) 
is feasible for a general case. For simplicity, so that 
[Ca^*]i(r) can be evaluated analytically, the excess buf- 
fer approximation is assumed (Smith, 1996) for our 
rough estimation, so 

[Ca2+ ]. (r) = [zc, / (zc,F2nDr) ] exp(-r /X), O) 

where X, the characteristic length, is {D/konlhY^'^. From 
parameters used in Wagner and Keizer (1994) , Smith et al. 
(1996), and Falcke (2003), we determined a high and 
low estimate for X. as 440 and 55 nm, respectively, and es- 
timated the range of [Ca^^]i(r) using Eq. 7 (Fig. 6). 

For r> 200 nm, [Ca^^]i(r) calculated from the two esti- 
mates for X. differs by over two orders of magnitude 
(Fig. 6), and it is no longer meaningful to consider 
[Ca"^]i(r) for a general case. More specific values for kon 
and -Bp are needed to give a better evaluation of [Ca^^Ji. 

From our estimation, it is clear that as [Ca^^]i(r) de- 
creaises for larger r, there is a range of r within which the 
Ca^^ released by an open InsPjR channel can activate a 
neighboring channel at that distance away, as long as 
[InsPs] is sufficiently elevated, for all InsPsRs studied 
(Foskett et al., 2007). Thus, the magnitude of ica ob- 
served confirms that CICR can be a mechanism to 
couple neighboring InsPgR channels to coordinate con- 
certed Ca^^ release by multiple channels to generate 
various intracellular Ca^^ signals. 

In this simple consideration, the [Ca^^Ji profile was 
estimated around an open channel with Po= 1 . In real- 
ity, the amount of Ca^^ moving through the releasing 
channel depends not only on «ca, but also on the sto- 
chastic gating of the releasing channel, which is dynam- 
ically regulated by the local [InsPj] and [Ca"*]i (Foskett 
et al., 2007) , which in turn are affected by Ca"^ released 
by the releasing channel and any neighboring active 
channels. Furthermore, the activation of an InsPjR 
channel by CICR is a complex dynamic process regu- 
lated by stochastic binding and unbinding of InsPj and 
Ca^^ to activating and inhibitory sites, which are affected 
not only by local [Ca"*] i, but also by the on- and off-rates 
of the sites that can be allosterically coupled (Atri et al., 
1993; Tang etal., 1996; Kaftan etal., 1997; Moraru etal., 
1999; Dawson et al., 2003; Mak et al., 2003) . To properly 
take into consideration all of these complicated factors 
affecting the regulation of Ca^^ signaling, quantitative 
kinetic modeling (De Young and Keizer, 1992; Swillens 
et al., 1994, 1998, 1999; Dupont and Swillens, 1996; 
Tang et al., 1996; Falcke et al., 2000; Sneyd and Falcke, 
2005; Swaminathan et al., 2009) using the right param- 
eters, of which ?ca is a critical one, is necessary. 

In summary, we have described the first direct electro- 
physiological measurements of the zca's driven by physi- 
ological [Ca^*] gradients across single InsPsR channels 



in a native ER membrane environment under physio- 
logical ionic conditions. These measurements will en- 
able more accurate evaluations of the amount of Ca^* 
released in a fundamental Ca"^ release event mediated 
by a single InsPsR channel, contribute to a better esti- 
mation of the coupling between the activities of neigh- 
boring InsPsR channels through CICR, and provide 
insights to improve future understanding and model- 
ing of intracellular Ca~* signals. 
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